Infrared self-consistent solutions of spinor QED 3 



o 

PL, 



Ph. 
43 



> 

(N 



X 



Tomasz RadozyckQ 
Faculty of Mathematics and Natural Sciences, College of Sciences, 
Cardinal Stefan Wyszynski University, Woycickiego 1/3, 01-938 Warsaw, Poland 

Three-dimensional electrodynamics in the spinor (i.e. two-component) version is considered. With 
the use of the so called Salam's vertex, the infinite hierarchy of Dyson- Schwinger equations is turned 
into a set of four self-consistent equations for four parameters describing the infrared behavior of 
fermion and boson propagators. It is shown numerically, that this set of equations has solutions, at 
least for certain values of gauge parameter. For weak coupling (i.e. for heavy fermions) the values 
of all these quantities are found analytically. In the case of massless bare fermion, masses of both 
particles are generated leading thereby to parity breaking. 
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I. INTRODUCTION 



Studying gauge field theories in low-dimensional space- 
time allowed for investigations both perturbative and 
nonperturbative aspects of quantum field theory. Among 
abelian models of this kind the particular mention de- 
serve quantum electrodynamics in two dimensions — the 
so called Schwinger Model (SM) — and in three di- 
mensions. The former, due to its full solvability, at least 
in the massless case, has become particularly fruitful and 
many nonperturbative effects have been discovered and 
analyzed. On can mention here the generation of the 
gauge boson mass, via the so called Schwinger mecha- 
nism, fermion confinement, chiral symmetry breaking, 
the presence of instantons and the ^-vacuum (see for 
instance @-[(|), bound state formation 0] but also the 
status of the perturbation series [8MTo| . Due to these 
nontrivial features Schwinger Model has become a test- 
ing laboratory for various aspects of the theory of strong 
interactions. The massive version of this model in no 
longer solvable but is still attractive thanks to its rela- 
tive simplicity and remarkable physical content [ll|> 03 ■ 

Quantum electrodynamics in three dimensions - the so 
called Planar QED — plays an important role in under- 
standing some aspects of QFT too. One can recall here 
confinement fl1i-[l7|]. chiral symmetry breaking (IH Il7l - 
[2pj . bound states j2l| - [23| and the analysis of the gauge 
dependence of the nonperturbative results with various 
approximations [24|,[25|. Contrary to the SM, the photon 
becomes in QED 3 a dynamic particle and corresponds to 
real degree of freedom, although it has some peculiar fea- 
tures connected with the fact that the magnetic field is 
now a (pseudo) scalar and electric one a two-component 
vector restricted by the Gauss law. The number of di- 
mensions attributes to this model also other specific fea- 
tures. The model may be formulated in two inequivalent 
versions: the fundamental fermion field may be chosen 
as four- or two-component one. In the latter — which 
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we concentrate on in the present work — the most inter- 
esting is the appearance of the topological photon mass 
term (2fjj . This term may be, just in three dimensions, 
added to the standard Lagrangian from the very begin- 
ning, since it exhibits both Lorentz and gauge (up to 
the surface terms) invariance — such a theory is called 
Chern-Simons QED3 — but even if absent, it is generated 
by the interactions with fermions. 

Both low-dimensional models find their applications 
also be yon d QFT: for instance in condense matter 
physics [13-[13- 

Many of the above mentioned phenomena, as con- 
finement, bound states or topological aspects, require 
for their study nonperturbative language. The prob- 
lems with perturbation calculations in QED3 manifest 
themselves also through ambiguity in proper regulariza- 
tion (26[. The set of Dyson- Schwinger (DS) equations 
- which, as one believes, contains the whole informa- 
tion about the quantum system — constitutes the main 
tool in nonperturbative investigations of QFT. The com- 
mon problem in this approach is, that these equations 
form an infinite hierarchy, the solution of which seems a 
hopeless task. To deal with this set one is forced to trun- 
cate this hierarchy at some level. In two dimensions (in 
SM), thanks to two gauge symmetries of the theory, it 
turned out to be possible to express the three-point (ver- 
tex) function by fermion propagator (i.e. two-point func- 
tion) and obtain the close equation for the latter. Not 
a whole hierarchy, but just one equation, which turned 
out to be easily solvable [3(| EJ. Similar method, but 
now applied to higher functions (five-point), allowed to 
write a self-consistent equation for two-fcrmion Green's 
function (32|. 

In greater number of dimensions, i.e. in QED3 and 
QED4, gauge symmetry is, however, not sufficient to 
fully express fermion-photon vertex function by lower 
ones and therefore we are obliged to assume a certain 
ansatz for T fJ, (p+k,p). Gauge covariance constitutes here 
a hint but does not solve the whole problem of truncat- 
ing the hierarchy. Various such ansatze have been pro- 
posed based on kinematic properties, renormalizability 
and requirement of satisfying the Ward-Takahashi iden- 
tity. One can mention in this context the Ball-Chiu [33j . 
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Burden-Roberts [34[ and Curtis-Pennington [35[ vertices 
with further improvements [36l [37j | . In the present work, 
however, we pay particular attention to the vertex intro- 
duced by Salam [38| and applied afterwards in the, so 
called, gauge technique (39l - l42| . The advantage of using 
Salam's ansatz in the form (|16j) . apart from its elegance, 
is that Feynman integrations in DS equations, are sim- 
ilar to those of perturbation calculations for which the 
computational technique is well elaborated. 

The main idea of the present work is similar to that of 
our previous paper 43] on QED 4 , where it was success- 
fully applied and the correct infrared form of propagators 
was established without the need of infinite renormaliza- 
tion. Here again we assume the possible infrared behavior 
of propagators up to a set of certain unknown constants 
and try to determine these constants from the require- 
ment of self-consistency. One has to emphasize that in 
this approach loop integrations are performed with in- 
frared forms of the Green's functions as integrands. Since 
we do not wish to extend the results beyond the infrared 
domain, this procedure is fully acceptable. Firstly we 
have to note that right hand sides of DS equations (fT8|) 
and (|28p have the form of convolutions, so the infrared be- 
havior of the integrals should be dictated by the infrared 
behavior of integrands, secondly our method should not 
be worse than ordinary perturbation calculations, where 
free functions are used as integrands. Finally one has to 
stress, that Salam's ansatz (I16|) in the infrared domain 
becomes exact. 

This work is organized as follows. In Section [TT] we 
collect the main definitions and properties of the model. 
In Section [IIII we formulate the assumptions as to the in- 
frared behavior of the photon propagator, fermion propa- 
gator and fermion-photon vertex. The following two sec- 
tions are devoted to obtaining, from the DS equations, 
the set of four nonlinear equations for four parameters: 
/3, r, Sm and Z3. In the last section we present the nu- 
merical results and the dependence of all the parameters 
on the strength of the interaction. We find also the ana- 
lytical results for the case of weak coupling. 

II. BASIC DEFINITIONS AND IDENTITIES 

Below we briefly summarize conventions used in this 
work and specify basic properties of Dirac gamma ma- 
trices in three dimensions. As to the set of matrices 7^ 
we use the two-dimensional, i.e. spinor representation in 
which all gamma's may be expressed by Pauli matrices 
in the following way: 

7 =*3= ^ o _,) , 7 =^=[i o ) , 
7 2 = i<?2 = (-1 ) ' W 

and matrix 7 s does not exist. Due to this nonexistence 
the phenomenon of chiral symmetry breaking can be in- 



vestigated only in four-dimensional representation, unless 
we introduce even number of flavors in the model. 

The metric tensor g pi/ and the totally antisymmetric 
tensor e^ a are defined as follows: 

.9 00 = -9 11 = -g 22 = 1 , 



with all other nonzero elements of the latter obtained by 
cyclic permutations. 

In three dimensions the peculiar feature of the 
choice ([1]) is that, apart from ordinary relations: 

{y*.7"} = *r, tr 7 ^ = 0, Tr [7^1 = 2^, 

Tr [yVW] = 2 (9^9 pa - 9 W 9 V ° + 9^9 vp ) , (3) 

also the trace of the product of three (i.e. odd number) 
gamma matrices has a nonzero value: 

Tr [7 M 7V] = -2ie^ p . (4) 

This identity is (mathematically) responsible for the pho- 
ton mass generation by the vacuum polarization loop. 

In sections IIVI and [V] we will also make use of other 
identities, which can be easily derived: 

e^T/iW = 2e^ Q 7 „ + 2i 7 V , 
e pva l^k a = -2i}t , (5) 

7^7„ = -#- 

Lagrangian density of three-dimensional electrody- 
namics with gauge fixing term has the following form 

C(x) = *(x) {vfd,, -m - eo^A^x)) V(x) (6) 
~F>»(x)F lu ,(x)-±(d lt A>>{x)) 2 , 

where A is the gauge parameter. The quantities mo and 
eo are the bare fermion mass and the bare coupling con- 
stant (charge) respectively and the former may eventu- 
ally be put equal to zero leading to the parity invariant 
theory (at least on the classical level). It is worth men- 
tioning that coupling constant has a dimension of v / mass 
and, as a result, the quantum theory is superrenormal- 
izable. The strength tensor F^ v has only three indepen- 
dent elements, two of which (F 10 and F 20 ) constitute 
to two-component vector of electric field, and one (F 12 ) 
pseudoscalar magnetic field. 

It is well known [26| . that in three space-time di- 
mensions to this Lagrangian density may be added 
the so called Chern-Simons (CS) term, proportional to 
e pva A^d v A a . Although it is not invariant under gauge 
transformations, its change reduces to a total derivative 
inessential for physical observables. Due to this term 
the photon field acquires a nonzero mass. In this work, 
however, we do not include explicitly CS term into the 
Lagrangian and the photon mass will be generated dy- 
namically. 
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III. ASSUMPTIONS FOR THE FUNDAMENTAL 
GREEN'S FUNCTIONS 

Underneath the low momenta forms of the dressed two- 
point Green's functions and three-point vertex are given. 
They reflect their analytic structures suggested by the 
perturbation calculation. 



A. Boson propagator 

For the free propagator of massless vector particle we 
have the usual formula: 



W 



WW 



1 WW 



A (fc 2 



(7) 



It is well known from the perturbative approach [26l . l44j , 
that in three space-time dimensions photon" field ac- 
quires a (topological) mass, due to the vacuum polariza- 
tion process. Even if CS term is initially absent in the 
Lagrangian, it is generated through the interaction with 
fermion loop. In the first approximation we can then 
assume for the dressed propagator 



D^ik) 
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WW 



(8) 



iZ 3 r 



k 2 (k 2 



2 ) 



z» vp k p - 



1 WW 



A (k 



2)2 



where r is the photon mass to be determined from the 
consistency conditions, with the renormalization con- 
stant Z 3 present only in the transverse part of D^ v {k). 
This is a known fact by virtue of the gauge invariance 
and the Word-Takahashi identity: 

k ll D^(k) = k ll D^"(k)=~ 1 ^, (9) 

which ascertains that only transverse part is modified by 
the interactions. 

By the reason of unitarity we expect the nonperturba- 
tive value of Z 3 to satisfy the condition < Z 3 < 1 [45j . 
In our work on QED4 we obtained this value to be equal 
to I Q. As we will see in section IVI1 the above bounds 
will also be satisfied in the present work. 

We treat the denominator k 2 —r 2 in (|8]) as the first term 
in its Taylor expansion around k 2 = r 2 and potentially 
admit also higher order polynomial in k 2 . 

The inverse of is 



D~ l (k) 



(10) 



(-* V + k^W) - iTe^ p k p - XWW 



B. Fermion propagator 

The free fermion propagator has the standard form: 

1 



- m 



(11) 



with a pole in the bare mass mo- Contrary to that, the 
perturbative calculations [26[ (at least in certain gauges) 
show that, even if initially massless ,, photon" acquires a 
mass r, the ,, dressed" fermion propagator should have, 
as a result of infrared divergences, a branch point at p 2 = 
to 2 , where to is a physical mass. Therefore we assume 
for further calculations the following form: 



S(P) 



1 



m)(l — p 2 /m 2 )P 



(12) 



with exponent (3 to be determined (together with mass 
renormalization constant 5m = m — mo) by the later re- 
quirement of consistency. The additional power 2/3 of p in 
denominator improves the high momentum behavior of 
loop integrals (after the analytical continuation of (fT2|)). 
making thereby the theory free of any ultraviolet diver- 



gences. 



Imz 





Rez 



FIG. 1: The contour of integration in formula (|13[) . which 
leads to the spectral function (1141) . 



To get the spectral representation of (fT2"j) . which will 
be needed for the construction of the vertex, we apply 
the Cauchy integral formula for the function 

/(*) 



(m 2 — z 2 )P ' 

choosing the contour as shown in Figure [T] deformed to 
sandwich both cuts on the real axis. This is possible for 
< P < 1, since otherwise the integrals over the small 
circles around branch points at z = ±to diverge. We 
have then 



1 



(m 2 — w 2 ) 13 



1 

c 

sin 7T/3 



dz 



(to 2 



z 2 )P{z- 
dx 



— m 

I 



[x 2 — m 2 )P(x — w) 



dx 

(x 2 — m 2 )fi{x — w) 



(13) 
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Applying this formula for w —p 1 and taking into account 
the additional factor to) in denominator of (fT2)) , we 
get that way (c.f. [H]) 



p(M) = 



sin(7r/3) 



7T (M -m)(M 2 /m 2 
x [9(M - m) - 9(-M - to)] , 



(14) 



where 9 is the Heaviside step function. We found it 
useful to separate out in the above definition of p(M) a 
term, that describes free propagation of fermion so that 
S(p) is represented as 



S(p) = / dMp(M) 



1 



1 



p 1 — to 



M 



(15) 



The equations, we obtain in the following sections, turn 
out to be self-consitent without assuming any nontrivial 
value of the fermion field renormalization constant Z 2 
and therefore, in the present work, we have just put it 
equal to unity. Eventual inclusion of this constant (i.e. an 
additional unknown variable) would require to generate 
more equations by adjusting higher Green's functions. 
This complicates our method, but we keep it in mind 
and leave for the next paper as a possible following step. 

The problem with properly defining constant Z 2 raised 
in [26j j does not arise in our approach since we do not 
require the full fermion propagator to have a pole at mass 



J 



to. On the contrary we know, that such a pole does not 
exist in p^|) due to the infrared divergences and only 
the additional coefficient in numerator come into play. 
Its value would be fixed by self-consistency and not by 
normalization requirement of the residue. 



C. Fermion-boson vertex 

Having expressed the full propagator S through the 
spectral function, we can now make use of the slightly 
modified (by inclusion of the second term) Salam's ansatz 
for the vertex function T^(p + k,p): 



S( P + k) T^(p + k,p)S(p) 
dMp(M) 

1 



• — TO p — TO 



+ ji-M 



r 



M 



(16) 



By acting on it with a vector k^, one automatically ob- 
tains the appropriate WT identity, which must be incor- 
porated in any method, that is expected to satisfy gauge 
invariance: 



k ll S(p + k)r' 1 ( P + k,p)S(p) = / dMp(M) 



= / dMp(M) 



1 



Ji- 



1 



1 



1 



+ jt-m j}-m p+ J6- M i>- M 
1111 



(17) 



+ )i-m i>-M ji + }t - M 
I 



= S(p)-S(p + k) . 



It is a necessary condition but does not guarantee by 
itself the invariance of physical quantities. 

The longitudinal part of the vertex is then correctly 
fixed, and transverse part (obviously (|T6| is not purely 
longitudinal) is postulated in an elegant way, naturally 
being a subject of further improvements (4ll. |42|| . 

The expressions ([l0]l. ([15]) and ([IB]), together with (fi"4]) 
will be now required to satisfy DS equations for small 
momenta. 



IV. DYSON-SCHWINGER EQUATION FOR 
THE BOSON PROPAGATOR 



Dyson-Schwinger equation for the photon propagator 
may be written as 



xTr 7o 



fia 



d 3 p 
(2tt)* 



k>*k a 1 k»k a 



k 2 A k 2 



S»-ie 2 x 



S(p)Tp(p,p-k)S( P -k)D^(k) 



(18) 



and has the graphical representation shown in Figure [2] 

W\AA/ = wwv + vwQm 

FIG. 2: Dyson-Schwinger equation for the gauge boson prop- 
agator _D M ^(fc). Light lines represent free propagators and 
heavy ones dressed propagators. The full circle stands for the 
full fermion-boson vertex. 



5 



This equation, rewritten for inverse propagator, will be 
more suitable for deriving self-consistency relations since 
D pv decouples from other functions. If we make use of 
the ansatz (fll)]) . it can be given the form 



D- 1 {kY v = -k 2 g pu + k»k v - Xk^k" + ie 2 Q Tr^ x 
X |,M„M,/|^(- > - ' 

1 ., 1 



m + ie 



7 



M + ie i> - M + ie 



m + ie 
(19) 



The expression under the integral is well known from 
one-loop perturbation calculation and its evaluation is 
straightforward, for instance by passing into euclidean 
space. It is needless to show any details, so we give only 
final (Minkowski) result. If we denote 



iC(*0 = 



(20) 



ielTijV 



d 3 p 



1 



1 



(27r) 3 p 1 — m + ie p 1 — ft ~ m + ie ' 
and introduce Feynman parameters, we find 



2tt 



(-JfeV + k"k u ) / dx 



x(l — x) 







' (m 2 - k 2 x(l - x)) 1 / 2 



1 r i 

--rnr^Pk / d-r 

2 "J (m 2 - fc 2 x(l - z)) 1 ^ ' 

o 



(21) 



where x is the Feynman parameter. When k 2 < Am 2 
both integrals are well defined. 

The tensor Hl^'(k) turns out to be transverse with 
respect to as is required by gauge invariance, and 
particularly by WT identity (j9|). What is important, is 
the appearance of the novel term, peculiar for three di- 
mensions, which is proportional to e pvp — a result of 
the nonzero trace of the product of three gamma matri- 
ces (H]). This property is responsible for the gauge field 
acquiring a mass. 

To find the full form of the polarization tensor we have 
to calculate — according to (fT6)l — the following integral 
involving the spectral function p(M): 



n^(fc) = / dMp(M) (ic(fc) - n£(*)) 



(22) 



Various integrals of that kind, containing spectral func- 
tion p(M), are collected in the appendix lAl When we 
apply formulae (|A3[) and (IA4[) . the only integrals that 



are left, are those over Feynman parameter x: 
e 2 r(/3 + l/2) ^ 



IP"(fc) 



27r 3 / 2 r(/? + l)m 







--mF^Pk f d-r - 

2 P J l-k 2 /m 2 x(l-x))P+V 2 

o 



(23) 



Both turn out to be expressed through the hypergeo- 
metric (Gauss) function (see appendix IBl formulae (|B1[) 
and (|B2[1 ). These results allow us to bring the whole 
Dyson-Schwinger equation ([T9|) to the form 



(24) 



D~ l {ky v = -\k»k v + {-k 2 g^ + k»k v ) x 



-ie^ p k 



gW + 1/2) 
127r 3 / 2 r(^ + l)m 
eW + 1/2) 

47r 3 / 2 r(/3 + 1) 



2 F 1 (2,/3 + l/2;5/2;fc 2 /4m 2 ) 
2 f 1 (l,/3 + l/2;3/2;fc 2 /4m 2 ) 



(25) 



and adjust the first two terms in the infrared domain, i.e. 
close to the photon mass (k 2 ~ r 2 ). In that way we get 
two equations for Z3 and r: 



Z3 1 = 1 



e 2 r(/3 + l/2) 



tZ. 



127T 3 / 2 77ir(/3 + 1) 

x 2 Fi(2, p + 1/2; 5/2; r 2 /4m 2 ) , (26) 
e 2 r(/3 + l/2) ^ 



3 47r 3 / 2 mr(^ + 1) 

x 2 Fi(l, p + 1/2; 3/2; r 2 /4m 2 ) , (27) 

which will be solved in section [VT] together with the other 
two derived form the DS equation for the fermion prop- 
agator. 



V. DYSON-SCHWINGER EQUATION FOR THE 
FERMION PROPAGATOR 



Dyson-Schwinger equation for the fermion propagator 
has the form 



S(p) 



1 



m 
d 3 k 



1 + ie 2 ol p x 



(28) 



■S{p + k)T v {p + k,v)S(p)D^{k) 



and may be represented as shown in Figure [3l 
It can be rewritten in a simpler manner 

(p'-m )S(p) = l + Z(p)S( P ) , 



(29) 
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_ . _ + r\ 

FIG. 3: Dyson-Schwinger equation for the fermion propaga- 
tor. As in Figure [2] heavy lines stand for full functions and 
light for free ones. 

where 

E(p)S(p) = (30) 

/d 3 k 
— S(p + k)T»(p + k,p)S(p)D,Ak) ■ 

I 



Now we have to substitute (fTBj) for the vertex, and (|SJ) 
for D^ v , obtaining 



E(p)S(p) = [E(p)£(p)U + [S(p)5(p)] B 

+[E(p)5(p)]o + [£(p)S(p)]z> , (31) 



where, to avoid lengthy expressions, we separated the 
contributions coming from different tensor structures in 
the photon propgator: g^" , k^k^, e fil/p k p , and from 
gauge-dependent longitudinal part: 



[£(p)S(p)U = 
-iZ 3 e§ y dMp(M) 

\E(p)S(p)] B = 

iZ 3 e 2 J dMp(M) 
p(p)5(p)] c = 



(2^r)3 



1 



1 



(2tt) 3 



+ ft — m + ie (k 2 — t 2 + is)(p' — m + ie) 
1 „ 1 



(m M) 



-Z 3 e 2 T J dMp(M) J 



d 3 k 



+ ft — m + is (k 2 + is)(k 2 — t 2 + ie)(jb — m + is) 

l k„ 



+ ft — to + ie (k 2 + is) (k 2 — t 2 + is)(p' — to + is) 
1 .. 1 



P(p)S(p)]d 

= 2 



d 3 k 
f 



- (to -> M) 

(m M) 



(32) 
(33) 
(34) 
(35) 



ft 



+ ft — m + is (k 2 + is) 2 (p' — m + is) 



(m -> M) 



Each of the above momentum integrals can be performed in an ordinary way by performing Wick's rotation and 
introducing Feynman parameters. In the symmetric integration no infinities arise even if the simple power counting 
indicate an apparent logarithmic divergence. Disregarding all details of the computation, as a result we have (again 
in Minkowski space) 



I A = ie 



I B = iel 



d 3 k 



1 



(27r) 3 p'+ft — m + ie tl k 2 — r 2 +is 8tt J (m 2 x — p 2 x(l — x) + r 2 (l — x)) 1 / 2 

o 

l 



f 



dx- 



p"(l — x) — 3m 



(36) 



d 3 k 



1 



(2tt) 3 * ft-m + is r {k 2 + is)(k 2 - t 2 + is) 8tt J ^[(m 2 x - p 2 x(l - x) + r 2 (l - x)) 1 / 2 



° / dx 



1 



Ic = e 2 Q 



+ to) 

V 2 
d 3 k 



(m 2 x — p 2 x(l — x)) 1 / 2 (m 2 x — p 2 x(l — x) + r 2 (l — a;)) 1 / 2 

l 



(p 1 — m) 



(37) 



(2tt) 3 

i>($ — to) 



l 



+ ft-m + is lu {k 2 +is)(k 2 - t 2 + is) Att 
x x 



dx 



-1 



T ■ £ 

Id = ie 



d 3 k 



(m 2 x — p 2 x(l — x) + r 2 (l — x)) 1 / 2 (ra 2 x — p 2 x(l — a:)) 1 / 2 

l 



(m 2 a; — p 2 x(l — x) + t 2 (1 — x)) 1 ' 2 

(38) 



ft- 



1 



1 



(27r) 3 p + ft - m + is {k 2 + is) 2 8n 
+ to) x(1-x) 

[p — m) 



° / dx 



1 



(m 2 x — p 2 x(l — x)) 1 / 2 



2 [m 2 x — p 2 x(l — x)) 3 / 2 



(39) 
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Now the contributions to T,(p)S(p) may be written as 



P(p)S(p)U = 

-Z 3 J dMp(M) 
\E(p)S(p)] B = 
Z 3 J dMp(M) Ij 

[£(p)S(p)]c - 
-Z 3 r y dMp(M) 

-j J dMp(M) 



Ia- 



1 



m + ie 



(m M) 



— m + 1£ 



■ — m + i£ 



(m -> M) 



(m -» M) 



(40) 



(41) 



(42) 



(43) 



1 — m + ze 



- (m -> M) 



and the appropriate spectral integrals, after a small rear- 
rangement, are carried out according to formulae of the 
appendix |XJ 

Before performing x integrations, it is important to 
notice, that if one substitutes ([12")) into the left hand side 
of (|28|) and multiply both sides by (tf+mo), there appear 
two singular terms, when p 2 — > to 2 : 



m 2 ^5m(p' + to) to 2/3 



(to 2 — p 2 )' 9 + 1 



(to 2 — p 2 )P 



(44) 



For our goal it is sufficient to pick out the identical terms 
from the right hand side of (|28|) . It is a very pleasant 
observation that it actually displays identical analytical 
structure, which is a strong sign of self-consistency. Con- 
sequently we do not need to preform explicitly all param- 
eter integrations, and can limit ourselves to those terms 
that contain singularities of the kind (|4"4")) . The appropri- 
ate results are given below: 



[£(p)S(p)U 



4Z; 



3 /C 1 (p 2 ,r 2 )x 



8tt 
+ to) 



(to 2 — p 2 )(l — p 2 /m 2 )P ' 
[E(p)S(p)] B « 

[mS(p)]c^-^flC 2 ( P 2 ,r 2 )x 



(45) 

(46) 
(47) 



(to 2 — p 2 )(l — p 2 /m 2 )P 



[S(p)S(p)] 



1 



87rAm \ (to 2 — p 2 )(l — p 2 /m 2 )@ 
1- ft 1 



2/3 to 2 (1- P 2 /to 2 )^ 



^(^ + to) 



(48) 



where w means that the expression does not diverge, 
when p 2 — > to 2 . The functions ICi and /C 2 have the fol- 



lowing form: 

* r Zl (49) 



1 

/C 2 (p 2 ,r 2 )- / dx 



( P 2 X 2 +T 2 (1- X)) 1 / 2 



(50) 



Equating identical divergent terms on both sides of 
Dyson-Schwinger equation (|29|) we obtain two relations 
for parameters Sm and ft: 



l = f|^f,S 2 (m,r)+ ' 



8tt 



AZ3to/3 



(51) 
(52) 



where auxiliary quantities S\(m, r) and 5 2 (m, r) are de- 
fined as 



Si=to/Ci(to 2 ,t 2 ) + 2t/C 2 (to 2 ,t 2 ) 



(53) 



«S 2 =/Ci(to 2 ,t 2 ) + 2to 2 



9/Ci(to 2 ,t 2 ) <9/C 2 (to 2 ,t 2 ) 



dm 2 



AmT- 



dm 2 



(54) 



and found in appendix [Bj 



VI. RESULTS AND CONCLUSIONS 

Below we rewrite the whole set of equations using 
renormalized quantities: fermion mass to, gauge coupling 
1/2 

constant e = Z 3 eo and gauge parameter A_r = Z3X. 
We additionally introduce a dimensionless parameter 
C= 3^5. obtaining 



<5to £ 



1 



2r 



x ln(2m/r + 1) + 



to 2to 2 ) 
1 ' 



i? 



( 


2H 


2r 




►4) 


2 




TO 


\ TO 


TO / 



x ln(2m/r + 1) + 



(55) 



(56) 



Z 3 = 



I - i T{ JL + } ,2 \ 2 *i (2, ft + 1/2; 5/2; t 2 /4to 2 ) , 



T 
TO 



c 



3 ^r(ft + 1) 
rCs + i/2) 



(57) 



-^-^ 2 F 1 (l,/3 + l/2;3/2;r7V). (58) 



Now the question arises whether this set of highly non- 
linear equations for unknown Sm /to, /3, and t/to has 
a certain domain of solvability. The answer is positive 
for a large range of values of gauge parameter Ai? and 
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for weak coupling. For exemplary value Xr — 1 numeri- 
cally found solutions of all equations as functions of the 
parameter £ are presented below. 

In the first graph of Figure U we show the dependence 
of the topological photon mass r (in units of to) on the 
parameter One can see, that for small values of pa- 
rameter it approaches the value e 2 /47r In spite of 
initial zero value of the photon mass, it reappears as an 
outcome of interactions with fermions. 

On the second graph the dependence of power f3 on £ is 
presented. As we remember our method is reliable only 
for < (3 < 1. This corresponds roughly to < £ < 1. 




FIG. 4: The dependence of photon mass r in units of m 
(upper plot) and power f3 (lower plot) on the parameter 
The dashed line corresponds to r/m = £. 

Figure [5] shows the fermion mass renormalization 5m 
(again in units of to) and coupling constant renormal- 
ization constant Z% as functions of £. Please note, that 
8m < to (for small £) and < Z3 < 1 and it decreases 
with increasing coupling. 

The comparison of all numerical results suggests that it 
might be possible to find analytical solutions in the case 
of weak coupling (which corresponds to heavy fermion) , 
i.e. when £ -C 1. From the Figure @] we conclude that 
in such a case we should assume r/m <C 1. Wc suppose 
also that /3 <C 1. All these assumptions will be justified 
a posteriori, when the whole set of parameters is deter- 
mined. 




FIG. 5: The dependence of mass renormalization 8m in units 
of m (upper plot) and charge renormalization constant Zz 
(lower plot) on the parameter 

For small values of r the equation (|56p may be given 
the approximate form (considering gauges, where £r is 
of order of unity) 



1 



C 1 



2 A fl /3 ' 



(59) 



since all other terms, particularly those containing a quo- 
tient t/to are small compared to 1//3 and have been ne- 
glected. Consequently we have 



P = 



2Xr 



(60) 



which is in agreement with the assumption /?« 1. 

For small values of r and {3 the equation reduces 
to (recall that T(l/2) = and 2 Fi(a, b; c; 0) = 1) 



TO 



(61) 



which is the Coleman's and Hill's result 44]. The first 
plot of Figure |4] shows, however, that for larger values 
of C it deviates form the straight line. This is well un- 
derstood since in Coleman's and Hill's paper the value 
of the photon mass was defined as r 2 = n(0), where n 
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is a polarization scalar. In our work the position of the 
photon mass is determined by the equation r 2 = II(t 2 ). 
One should also keep in mind that our solution is only 
approximate. 

From (|61[) we find again that the assumption r/m <C 1 
is justified. In a similar manner we get the following 
result for Z 3 : 



< Z 3 = 1 - j- < 1 



and for 5m: 



5m 
m 



l-21n(C/2) 



(62) 



(63) 



The unpleasant feature is the gauge dependence of the 
mass renormalization, which might indicate the similar 
gauge dependence of the physical fermion mass. It is, 
however, a common feature of perturbative and nonper- 
turbative results in QED 3 [H, |47j • 

In the weak coupling regime one can obtain also the an- 
alytical results for the Landau gauge, which corresponds 
to taking the limit Xr — > 00. As is known, in this case 
the infrared singularities disappear (2(| |4(| so we expect 
f3 — s- 0, since the fermion propagator should have an or- 
dinary pole at physical mass and not a branch point. 
Therefore, before taking the limit, we put /? = Po/Xr 
with constant Po to be determined. With all other as- 
sumptions identical as above, we get similar results as 
before: r/m = £, Z 3 = 1 — £/3 and 



-IK 



and hence for ( <C 1 



P 



c 



2A fl (l - C) 



Po 



2X R 



c 



2(1-0 



Xr — >oc 



(64) 



(65) 



The result for the mass renormalization in this case is 



5m 
rn 



c 



(l-21n(C/2)) 



(66) 



The nonanalyticity in £ observed in the last formula 
and in ("S3")) , which manifests itself through the presence of 
a logarithmic function of the coupling constant, is similar 
in nature to the one found by us in QED4 [43J, where it 
constituted a reflection of the possible ill behavior of the 
perturbation series [Hj]. 

The solutions of the set of nonlinear equations ([55] - [58j) 
involving gauge parameter Xr can, in general, be gauge 
dependent. To clarify this point we performed a couple 
of plots of the model parameters for different gauges, i.e. 
for A^ = 1,3, 10,oo, the last case corresponding to the 
Landau gauge. In Figure [5] the dependence of r and j3 on 
Xr is demonstrated. Figure[7]shows the same dependence 
of 5m and Z 3 . 

It may be easily seen, that, apart from the parame- 
ter ft, which is and may be obviously gauge dependent, 



0.5 




FIG. 6: The comparison of the behavior of photon mass r 
(upper plot) and parameter /3 (lower plot) for different gauges: 
solid line - Xr = 1, dotted line - Xr — 3, dashed line - 
Xr = 10, mixed line - Xr = 00 (Landau gauge). The latter is 
not visible on the upper plot, since it is overlapping with the 
solid line. 



all other values become practically gauge independent 
for weak coupling. To get the reliable results for strong 
coupling we need to go beyond the first approximation 
consisting of the simple assumptions ©, (ll"2"|) and (fT6"l) . 
The same refers to other nonperturbative (but not exact) 
calculations. 

At the end we would like to note, that from the de- 
pendence of 5m (Figure [S] and Figure [7]) one sees, that 
for certain values of parameter £ we have 5m = m (i.e. 
mo = 0). Assuming value 5m /m = 1 as fixed, the equa- 
tions ([55)1 - (15*5)) may be in turn solved for (, t, /3 and Z 3 . 
In the case of the Landau gauge the following values are 
obtained 



r 
rn 



C « 0.34 



(67) 



from which we get the generated masses (the Lagrangian 
in this case does not contain any masses) to be 



2.94 x — , 

4-7T 



e 



(68) 



and the parity becomes broken. For two other parameters 
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we find 

(3 « 0.21 (naturally f3 = 0) , Z 3 w 0.89 . (69) 

A comment should be made here. The obtained result 
of parity breaking is not at variant with Vafa-Witten the- 
orem [43, H3] stating, that there is no spontaneous parity 
breaking in theories with vectorlike fermions. The three 
-dimensional QED with two-component fermions is very 
special theory also in that aspect, that it evades the proof 
of this theorem. The crucial point in the Vafa's and Wit- 
ten's argument is the positivity of the boson measure 
(after integrating out the fermion degrees of freedom) 
which demands the positivity of the determinant of the 
Dirac operator. This positivity may be easily proved in 
four component version of the theory, because by the ap- 
plication of the matrix 7 s one sees that eigenvalues are 
always paired in such a way that the unwanted signs can- 
cel. However, in spinor version of QED no such 7 15 matrix 
exists (there is no matrix anticommuting with matrices 
7 , 7 1 and 7 2 ) so the above argument fails. 

The results of [5l[ , which support the Vafa-Witten con- 
clusion, are not in conflict with the present work either. 
They were obtained with the assumptions to <C a or 
to ^> a, which correspond to ( > 1 or ( < 1 respec- 
tively, while in our case we have £ ~ 0.34. As authors 
write: 'We have not yet excluded the possibility that a 
mass of order a is generated dynamically'. On has to 
stress also, that the conclusions of [5l| are got in the 
limit of infinitely many flavors, while we have one flavor 
theory. The number of flavors, particularly if it is even, 
may be essential for the positivity of the determinant 
spoken above. 




FIG. 7: The comparison of the behavior of fermion mass 
renorrnalization Sm (upper plot) and charge renormalization 
constant Z3 (lower plot) for different gauges: solid line - 
Xr — 1, dotted line - Xr = 3, dashed line - \r = 10, mixed 
line - Xr = 00 (Landau gauge). On the upper plot the lines 
for Xr = 3 and for Landau gauge follow close curves. 



Appendix A: Spectral integrals 

Below we collect the set of integrals involving spectral 
density (fl"4"|) . that we used in sections HVl and IVl First one 
can find that for a 2 < to 2 , one has 



/ 



dMp(M) 



( m 2_ a 2)l/2 ( M 2_ a 2)l/2 

m 2 ^r(/3 + 1/2) 



v /HX ( S + l)(m 2 -ay+ 1 /2 



(Al) 



and similarly 



dMp(M) 



(to 2 - a 2 ) 3 / 2 (M 2 - a 2 ) 3 / 2 
2TO 2/3 r(/3 + 3/2) 



v/HX/? + l)(m 2 -a 2 )£+ 3 / 2 



(A2) 



where T is the Eulcr function. From the former result by 
the appropriate substitution we obtain 



dMp(M) ( {m2 _ k2 l x{l _ x))1/2 - (m - M)) 

r(/3 + i/2) 



,2J 



^T{j3 + 1) (m 2 - k 2 x{l - x))^ 1 ' 2 ' 

(Tfl 

r(/? + i/2) 



(A3) 



TO 



2/3+1 



^T{P + 1) (to 2 - k 2 x{\ - x))P +1 / 2 ' 
where we additionally used the relation 



(A4) 



dM p(M)M f(M 2 ) = m J dMp(M)f(M 2 ) . (A5) 
With a slight modification of (|A3[) and (|A4[) one can 
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easily get the other set of integrals: 
dMp(M) ' 



(m M) 



(to 2 x — k 2 x(l — x)) 1 / 2 

r(/3+l/2) m^x 13 

v^r(/3 + 1) (m 2 x - k 2 x{l - x))^ 1 / 2 ' 



(A6) 



r(/? + i/2) 



(m M) 



0Fr(/3 + 1) (m 2 a; - fc 2 x(l - x))^ 1 / 2 ' 
Prom (|A2[) one finds 

1 



(A7) 



dMp(M) { {m , x _ k2x{1 _ x))3/3 ~ (m - M) 
_2r(/3 + 3/2) m 2 ^ 



0Fr(/3 + 1) (m 2 x - k 2 x(l - x))P+ 3 / 2 



(A8) 



^ ^ \ (m 2 x — k 2 x(l — a;)) 3 / 2 
2T(J3 + 3/2) 



0Fr(/3 + 1) (m 2 x - fc 2 x{l - x)Y+z/ 



(m M) 

2 ' ( A9 ) 



and also 

dMp(M) x 



(m 2 x - fc 2 a;(l - x) + r 2 (l - x)) 1 / 2 
r(/3 + l/2) 



(m M) 



m 2 Px l3 



V^r(/3 + 1) {m 2 x - k 2 x(l - x) + r 2 (l - x))P+V 2 ' 

(A10) 

The last integral, we will need, is lengthy and therefore 
we will not give it in its full complexity. Happily in Sec- 
tion |V] only singular terms when p 2 m 2 are necesssary: 



Appendix B: Parametric integrals 

The parametric integrals needed for the vacuum polar- 
ization tensor in Section IIVI are 



dx 



x(l — x) 



o (l-x{l~x)k 2 /m 2 )l 3 + 1 / 2 

1 F 1 (2,/3 + l/2;5/2;fc 2 /4m 2 ) . (Bl) 



G 



dx 



1 



o (l-x(l-x)k 2 /m 2 y+ 1 / 2 

= 2 F 1 (l,/3 + l/2;3/2;fc74m 2 ) . (B2) 



In Section [V] there appear two other parametric inte- 
grals, we denoted by Si and 52 (formulae ([55)) and ([54")) ). 
which may, after elementary calculation, be expressed by 
simple functions: 



Si = mJCi (to 2 , t 2 ) + 2t/C 2 (to 2 , r 2 ) 
l 

m(x + 2) + 2r 



(B3) 



= dx 



(m 2 x 2 + r 2 (l — a;)) 1 / 2 

2 





t ( 2r 
1 + I 2 + — 

TO 



, ,ln(2TO/r+l) 
m 2m 2 ' 



S 2 =/Ci(to 2 ,t 2 ) + 2to : 



,<9/Ci(to 2 ,t 2 ) 
dm 2 



ArriT 



g>/C 2 (to 2 ,t 2 ) 
dm 2 



t l dx 



(2m + t)x 2 + t(x - 2) 
[m 2 x 2 +T 2 {l-x)y/ 2 



1 




2r 




^) 


TO 


24 


TO 


\ TO 


TO / 



(B4) 



where K-i and /C 2 were defined in equations ([49)) and ([50 



rfMp(M) l^ X (AU) 



^ - (to ^ M) 



(to 2 x — j* 2 .t(1 — cc) + t 2 (1 — x)) 1 / 

^ TO 2 ' 3 1 

: (p 2 X 2 + T 2 (l - X)) 1 / 2 (TO 2 - p 2 )^ 1 
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